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FOREWORD 


This final report has been prepared in accordance with require- 
ments of Contract NAS5-11445 to present data and conclusions from 
a nine-month study for Goddard Space Flight Center by the Martin 
Marietta Corporation, Denver Division, The work was done under 
the management of the NASA Project Manager, Mr, George M. Levin, 
Advanced Plans Staff, NASA-Goddard Space Flight Center, The 
report is divided into the following volumes: 

I - SUMMARY 

II - SUPPORTING TECHNICAL STUDIES 
III - APPENDIXES 
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APPENDIX A 

HEAT SINK THERMAL STRESS ANALYSIS 
B. J. Buell & H. Hotchkiss 


November 10, 1971 



In a hemispherical probe shell entering Jupiter f s atmosphere at 
high velocities, the length of survival time is critical for ade- 
quate data recovery. A structural shell of beryllium is being 
considered for use and has been analyzed for thermal stresses at 
two altitudes during entry: 70 and 80 kilometers below turbo- 

pause. 

The best available method of analysis is the Wilson Axisymmetric 
Analysis (Martin Marietta Computer Program File Number SA005) . 

The program permits use of multiple bilinear stress-strain curves 
for various temperatures of the structural material and calculates 
stresses and deflections resulting from a temperature variation 
through the thickness, and ^ variation around a shell of revolu- 
tion by means of an iterative procedure converging on stress- 
strain pairs falling on or near the appropriate interpolated 
stress -st rain curve. 

The results of the -70 and -80 kilometer analyses are attached. 

The analytical model used (Figures A-l and A-2) consists of ele- 
ments and nodal points with corresponding temperatures. The 
beryllium material stress-strain curves and the respective bi- 
linear stress-strain curves used are shown in Figure A-3. The 
tangential stresses versus location calculated by the program 
are shown in Figures A-4 and A-5 for the -70 and -80 kilometer 
altitudes, respectively. The resulting calculated deflections 
for the two altitudes are presented in Table A-I. 

The maximum tensile stress in the inner row of elements occurs 
around 47° from the apex for the -80 km model and around 33° for 
the -70 km model. The calculated stress is considerably in ex- 
cess of the stress from the bilinear curve at the corresponding 
temperature, because the program calculates effective stress 
from an effective modulus of elasticity. Since the calculated 
maximum tensile strain (0.0014 in. /in.) is of reasonable magni- 
tude and the un.‘ axial stress-strain curves at elevated tempera- 
tures show the beryllium unable to sustain the calculated high 
compressive stress, the resulting high calculated tensile stresses 
are not expected to cause failure unless an unexpectedly severe 
reduction in ductility and strength should accompany the biaxial 
state of tensile stress. Multiaxial tensile stress-strain data 
should be accumulated or actual prototype testing of the beryllium 
sphere should be planned if beryllium is considered for use. 
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Figure A-l Elements and Nodal Points Temperature Distribution 
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Figure A-2 Elements and Nodal Points Temperature Distribution j 
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Figure A-3 Stress-Strain Curves for Hot-Pressed Beryllium 
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Figure A-4 Tangential Stress , for Outer and Inner Elements; 
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APPENDIX B 


SOLID PROPELLANT ROCKET MOTOR 
EXHAUST PLUME STUDY 


A. Reichert 


February 25, 1972 


The effect of solid particulate exhaust matter from the solid 
rocket deflection motor impacting the spacecraft could pose a 
problem for the spacecraft surface or science instruments. It 
was found that most existing exhaust plume computer programs 
were applicable for gaseous exhaust products only, and limited 
funding precluded the adaptation of a Martin Marietta computer 
program for solid rocket motor exhaust plume analysis. 

The solid particulate matter in a solid rocket motor exhaust 
plume remains within the projected cone of the motor nozzle, and 
when fired in a vacuum, maintains its nozzle exit velocity. To 
prevent this particulate matter from impinging upon the space- 
craft, a dual nozzle motor design was proposed; this design di- 
rects the particulate matter away from the probe axis which is 
aligned with the spacecraft. 

An alternative approach investigated was the ignition of a single 
nozzle (axial mounted) deflection motor. In this approach, the 
separation distance between the probe and spacecraft vas in- 
creased before deflection motor ignition. This method reduced 
the number of particles impacting the spacecraft as the distance 
increases, but does not remove the spacecraft from the path of 
the particles. 

Relative costs of single and dual nozzle motors and tolerance of 
spacecraft and instruments to small quantities of exhaust product 
impingement should be considered in greater detail before final 
design selection. 

Several solid propellant rocket motor manufacturers were contacted 
with respect to solid propellant exhaust plume characteristics in 
a vacuum. The concensus stated that the solid particulate matter 
would lie just within the rocket exhaust cone angle. This angle 
is generally equivalent to a 30° cone. A typical impingement 
analysis for the Turbopause Probe and Pioneer Spacecraft is pre- 
sented in section A; the difference in the cost quoted by three 
of the solid motor manufacturers contacted are presented in sec- 
tion B; selection of the solid motor design and reason for its 
selection is presented in section C. , 
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TYPICAL EXHAUST PLUME ANALYSIS 



1) Ai ■ - » 225,200 ft 2 (cross sectional, area of plume 

1000 ft aft of nozzle exit). 

2) Propellant weight in the motor is 6.09 lbm of which 30% 

is solid particulate matter; therefore 1.827 lb will be spread 
over 225,200 ft 2 * (or 81 x 10 -7 lbm/ft 2 ). 

3) Projected area of the spacecraft is assumed to be 100 ft 2 ; 
then the total propellant impacting the spacecraft is: 

81 x 10" 5 lbm 


(81 x 10” 5 ) (9 x 10 3 ) 

mi vi - m 2 v 2 , A Ufa L - v 2 

v 2 ■ 0.0133 ft/sec added to the spacecraft 
where 


mi * mass of particulate 
vi ■ velocity of particulate 
m 2 * mass of spacecraft 

v 2 • velocity impacted to spacecraft from impacting particulate. 




The average diameter of A&2O3 molecule — 

D ■ 3.16 x (0.337) 2 microns 
ave 

■ 0.358 microns “ 0.358 x 10~ 6 meters 
- 1.175 x 10“ 6 ft. 

The volume of a molecule ■ ~ (0.5875 x 10 6 ft) 3 ■ 0.85 x 10~ 18 
ft 3 3 

The mass of a molecule *» 247.6 x 0,85 x 10 18 - 2.1 x 10~ 16 lbm 

per AJI2O3 molecule. 

Based on the large proportion of AJI2O3 forming the solid particu- 
late matter of the rocket exhaust, assume the 81 x 10~ 5 lbm inter- 
cepting the spacecraft is AJI2O3. The number of molecules of 
AH2O3 then becomes: 

81 x 10 5 lbm molecule „„ , , n] 1 

38.6 x 10 11 molecules 

2.1 x 10" 16 lbm 

that will hit the spacecraft, or 


38.6 x 10 11 molecules 
100 ft* 


38.6 x 10 9 hits/ft 2 


The penetration of Mylar film by these solid motor exhaust parti- 
cles can be calculated as follows: 1 

c i p * 5 

where 


p ■ particle density (g/cc) 

T ■ target thickness penetrated (m) 
d ■ particle diameter (m) 

Ci ■ penetration coefficient (cc/g) 


B-3 


then T - (0.3) (2.7 g/cc) (0.36 x 10~ 6 meters) 
■ (0.81) (0.36) (10 6 meters) 

* 0.29 x 10" G meters 
or T ■ 0.29 x 10 3 mm 


Typical exhaust products that will be Included in those impact- 
ing the spacecraft are tabulated. 


Product quantity 
Exhaust product moles /100 gm 


SiO 0.0001 
A£C0 2 0.0006 
MC0 3 0.0001 
A£ 2 0 3 0.3334 


COST ESTIMATE COMPARISON OF OFF-THE-SHELF SINGLE AND DUAL NOZZLE 
SOLID ROCKET MOTORS (BY COMPANY) 


Hercules 


No off-the-shelf-motor. Cost of one inert (tooling) and two live 
motors delivered, plus development, is $125,000 for either single 
or dual nozzle design. 

Thiokol 


Single nozzle (TE-M-541) design off-the-shelf costs $3,000 to 
$3,500 per motor. 

Modification of TE-M-541 to dual nozzle configuration $50,000, 
plus test of 5 motors for $40,000 more. Delivery of motors for 
$3,000 to $3,500 each. 

UTC 

No off-the-shelf hardware. No difference in cost between dual and 
single nozzle designs. Cost not available. 

Mi arome te avoid I'ene ivcticm Test Data for Plastic Materials. 
TN-3770092 (Viking Engineering Report) , Martin Marietta Corpora- 
tion, February 12, 1971 
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SELECTION OF MOTOR 


A dual nozzle solid propellant rocket motor design was selected 
for application to the Jupiter Turbopause Probe for deflection 
velocity attainment. Such a design was selected because It as- 
sures there will be no solid particulate matter exhaust products 
Impinging on the spacecraft from the probe deflection motor. It 
Is also doubtful that an off-the-shelf motor (either single or 
dual nozzle) will be available for direct application to the 
probe without some modification. The cost of developing a single 
or dual nozzle motor for this application is the same (approxi- 
mately $125,000). The dual nozzles centerlines should be sepa- 
rated by an angle of approximately 45 degrees. 
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APPENDIX C 


PARAMETRIC STUDY OF 
DEFLECTION MOTOR PROPULSION SYSTEM 


A. Reichert 


* 


October 18, 


1971 


Parametric studies were conducted on three different propulsion 
systems covering a total impulse range from 500 to 25,000 Ibf-sec. 
The three propulsion systems studied were a single thruster pres- 
sure regulated GN 2 system; a single thruster hydrazine direct 
catalytic (liquid feed) system using a GN 2 blowndown pressuriza- 
tion system; and a single solid propellant motor system. All 
three propulsion systems delivered the required total impulse in 
a single firing. Propellant acquisition for the hydrazine system 
is not a problem as the probe is spin-stabilized prior to initiat- 
ing propulsion system thrusting. 

1. GN? Pressure Regulated Cold Gas Propulsion System 

The GN 2 propulsion system studied is shown schematically in Figure 
C-l and individual weights of the components comprising the feed 
system are indicated. The tank is initially filled to a pressure 
of 3500 psia, which will blow down to 200 psia during system 
usage resulting in a propellant residual of 5.7% by weight. 

Based on past performance on many programs, titanium was selected 
as the tankage material. Figure B-6 of Ref 1 presents a curve of 
mass of propellant and tank as a function of mass of propellant 
for a range of propellant mass from 0.2 to 10 4 lbm. Adding the 
weights of the components shown on Figure B-l to the tank plus 
propellant weight for a given propellant from the Ref 1 (Fig. B-6) , 
gives the initial propulsion system mass as a function of propel- 
lant mass. For steady state operation, the system will have a 
delivered of 72 Ibf-sec/lbm. The product of propellant mass 

and I gives I (total impulse). Figure C-2 presents the re- 
s p x 

suiting as a function of W pg (total propulsion system weight) 
for a range of from 100 to 2500 lbf-sec, 

2 . Single Thruster Hydrazine Direct Catalytic (liquid Feed) System 
Using a GN2 Blowdown Pressurization Feed System 

This system is shown schematically on Figure C-3. Included are 
the system component weights (less tank, propellant and pressurant) . 
Figure B-8 of Ref 1 presents a GN 2 blowdown feed system mass, in- 
cluding propellant and pressurant, as a function of propellant 
mass for a range of propellant mass from 0.04 to 400 lbm. Figure 
3-9 of Ref 1 presents a schematic of the feed system and component 
weights that are included in the feed system mass of Figure B-8 
of Ref 1. The difference in feed system component mass between 
Figure B-9 of Ref 1 and Figure C-3 was added to the feed system 
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Figure C-3 Deflection Propulsion System Catalytic Hydrazine (Blowdoun) 




mass of: Figure B-8 of Ref 1 to make the blowdown system curve 
apply to the Figure C-3 configuration. Figure C-4 presents the 
propellant fraction (X) of the Figure C-3 configuration as a func- 
tion of for a range of I from 500 to 30,000 lbf-sec. For 

this system, the was based on an average of 225 lbf-sec/ 

lbm of propellant. A survey of this type of blowdown system in- 
dicates that a will not exceed 2 lbf-sec/lbm. Selecting a 

constant of 225 lbf-sec/lbm will be a very close, average 1^ 

approximation. Figure C-5 presents the initial weight of the 
propulsion system as a function of (total impulse). 

Thrust will vary with thrust chamber inlet pressure as shown on 
Figure C-6. This effect must be considered when determining the 
time required to attain a given impulsive velocity increment, or 
magnitude of acceleration associated with the thrusting period. 
Figure C-7 presents tank top pressure as a function of percentage 
of propellant expended for initial tank cop pressures of 300 and 
400 psia. The tanks for the blowdown pressurization method are 
assumed to have an initial volume equal to twice that required 
for the propellant alone. The GN 2 pressurant is assumed to ex- 
pand in a polytropic process with the polytropic exponent equal 
to 1.29 (obtained from Titan IIIC Transtage ACS data). A pressure 
drop of 20 psi is assumed from tank top pressure to thruster in- 
let pressure. 
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Figure C-6 Thrust vs Inlet Pressure N 2 H^ Direct Catalyst Thruster 



(Bjsd) ajnseaaa dox ^fu»x 


*d 

•3 

8 

ft 


a 

o 

M 

eu 


C-9 


Figure C-7 Tank Pressure as a Function of Propellant Expended (Blmdom.) 



Single Solid Propellant Rocket Motor 


Solid rocket motor manufacturers were contacted to obtain recom- 
mended techniques for parametrically sizing solid rocket motors 
over a range of 1 T from 500 to 25,000 lbf-sec. Computer programs 

for rough sizing and weight determinations were available; how- 
ever they included proprietary information and therefore could 
not be furnished. 


It was therefore decided that the reference probe deflection 
motor propellant characteristics and motor geometric relation- 
ships would form the basis for a parametric study of solid pro- 
pellant rocket motors over the range of I T noted. Figure C-8 

presents the propellant and motor characteristics for a 1^ of 

1746 lbf-sec and Figure C-9 shows the geometric parameters studied 
as a function of 1^ of the motor. For purposes of this study, the 

thrust coefficient and chamber pressure were assumed constant as 

was the delivered I of the motor. As a result, the throat area 

sp 

varied directly as the thrust (F) varied, which in turn was 

directly proportional to the propellant weight / W_ ) and total im- 
pulse (Ij. ' P/ 


Figure C-10 presents the propellant fraction (X) as a function of 
I T< This figure was constructed from representative solid motor 

data. Figure C-10 is used to relate W as a function of I as 

Fb I 


presented in Figure Oil* Historically, it has been possible to 
achieve a higher propellant loading as the of the motor in- 
creases, as is indicated in Figure 012. This curve will tend 
to become asymptotic around propellant loadings of 0.85. The 
data presented is used in the calculation of motor dimensions as 
a function of 1^. Figure 013 presents W p ~ as a function of 1^ 


and is based on a I of 287 lbf-sec/lbm. Figure 014 presents 

F as a function of and is based on the assumptions of propel- 


lant and motor characteristics previously discussed. Figure 015 
presents chamber length as a function of 1^ and is based on W p as 

a function of I T> Vp^V^ as function of 1^, the geometric relation- 
ships of the reference vehicle motor, and the Figure 09 schematic 
The motor chamber diameter nozzle length ^L^j, and throat 

area (A.pjas functions of I T (Figures 016, 017, and 018, re- 
spectively) are based on the same data sources as was used for the 
chamber length. 


010 


Perf orman ce 

Total Impulse, lbf-sec 1746 

Average Thrust, lbf 318 

Average Pressure, psia *30 

Delivered I , lbf-sec /Ibm 4.66*5 

sp 

Duration, sec 5,48 

Propellant Density, lb/in. 3 0.063 

Reliability 0*998 


Weight Summary 

Total Inerts, lbm 3.58 

Propellant, lbm 6 .09 

Total Motor, lbm ' 9,67 

Mass Fraction 0.630 


Moment of Inertia 
Roll 40 lb-in. 2 

Pitch or Yaw 84 lb-ir. 2 



Figure C-8 Probe Deflection Motor 
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Figure C-9 Solid Propellant Rooket Motor Geometry 



C-13 





Total Inpulse, I T (lbf-sec) 


Figure C-ll I T vs foi* State -of- the -Art Solid Propellant Motors 


C-14 


Total Impulse, I T (lbf-sec) 



0.6 0.7 0.8 

Propellant Loading, Vp^V^ 

Figure C-12 Solid Propellant Motor Volumetric Solids Loading Fraction 
as a Function of Total Impulse 
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Figure C-16 Solid Rocket Deflection Motor D p vs 
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Figure 019 presents a comparison of the three types of deflection 
propulsion systems studied using W p g as a function of T^. As can 

he seen, total Impulse Is not the sole consideration that prompted 
selection of a solid rocket mcior for the deflection maneuver ve- 
locity increment. The exceptions could be associated with the 
hisrh thrust levels and single thrusi: applications of the solid 
motor. 


Reference? 


1. Lee B. Holcomb: Satellite Auxiliary - Propulsion Selection 

Techniques. NASA, JPL TR 32-1505. 
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APPENDIX D 


MONTE CARLO DEFLECTION DISPERSION ANALYSIS 


E. D. Vogt 


December 1971 


A Monte Carlo computer program is used to compute the dispersions 
In communication and entry parameters caused by errors and uncer- 
tainties at the time of the deflection maneuver. The deflection 
maneuver itself is defined in detail in subsection 1V.D.1 of 
Volume II of this report. A summary of the Monte Carlo technique 
is supplied in this appendix. 



A. 


ERROR SOURCES 


Two distinct types of errors are identified as causing dispersions 
from the nominal entry parameters. First, because of errors in 
the guidance and navigation process prior to the deflection man- 
euver, there are uncertainties in the spacecraft state at de- 
flection. Secondly, there will be execution errors made by the 
spacecraft and probe in the implementation of the required man- 
euver. 

1. Guidance and Navigation Uncertainties 

Quantitative measures of the uncertainties in the spacecraft 

state at the deflection point are provided by the control and 

knowledge covariances* at that point. The control covariance 

P is a 6 < 6 matrix defined by 
c 

P c ■ E [ ( X a c t - *h«) ( X .ct - X :*om) T ] 

where E is the expectation operator and X . is the random vari- 

act 

able vector describing the actual state (6 vector of position and 

velocity) of the spacecraft and X is the nominal state of the 

nom 

spacecraft. P £ thus gives a measure of tne probabilistic devia- 
tion of the actual deflection state from the nominal state. 

The generation of the control covariance proceeds as follows. 

It is assumed that the control errors result solely from errors 
at the last midcourse correction prior to deflection. A further 
assumption is made that the last midcourse correction is small 
enough that the execution errors are dominated by the knowledge 
errors at the time of the correction. A large a priori knowledge 
covariance is assumed at 25 days prior to the last midcourse 
maneuver when the tracking for that midcourse is initiated. The 
a priori knowledge covariance is reduced by processing simulated 
measurements for 25 days to determine the knowledge covariance 
at the midcourse. This covariance is then simply propagated 
(processing no measurements) over the last 13 days to generate 
the control covariance at the deflection point. 

*The control and knowledge covariances referred to in this report 
correspond to the correlation matrix X of the actual deviation 
vector and the correlation rating E of the estimated errors in 
Battin's notation of Reference 1. General analytic details of 
the formulation of these mathematical tools may be found in this 
reference. 
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The knowledge covariance at any point is defined mathematically 
as : 


P . - E f/x - 

k |^\ est 

where X ^ is the 
est 

states if X is 
act 

tative measure of 
of deflection. 



random variable representing possible estimated 
the actual state. Thus, P u provides a quanti- 


the estimation error to be expected at the time 


Tne knowledge covariance is generated similarly to the control 
covariance. Tracking begins at 38 days prior to deflection and 
measurement processing continues to one day prior to deflection. 
The knowledge covariance at that time (the deflection maneuver 
computation time) is then simply propagated to the nominal de- 
flection time. 


Thus, the knowledge and control uncertainties are really functions 
of the approach orbit determination (OD) processes. Since the 
orbit determination process is highly dependent upon the detection 
of the gravitational effects of Jupiter, the knowledge and con- 
trol uncertainties decrease as the deflection maneuver is per- 
formed nearer Jupiter. 

Figure D-l presents a summary of the control and knowledge covar- 
iance computations. The pertinent data used in generating the 
uncertainties is supplied in Figure D-1C along with a pictorial 
representation of the process. Deflection radii of 10- , 30- , and 
50-million kilometers were studied. These radii correspond to 
time intervals of approximately 8, 25, and 44 days before encoun- 
ter. For any of the deflection radii, the tracking was initiated 
38 days before deflection tracking through 25 days to generate 
the midcourse knowledge covariance. That covariance was then 
propagated forward to deflection to generate the deflection con- 
trol co^-riance. Tracking was reinitiated on the midcourse 
knowle ^ covariance and carried through to one day prior to 
deflection. This was then simply propagated over the final day 
to generate the deflection knowledge covariance. 

The dominant part of the deflection covariance matrice is the 
upper left 3><3 partition defining position uncertainties. 

This may be rotated into the standard RST coordinate system. 

Then the one sigma uncertainty in the spatial miss or the impact 
parameter B is defined in the TR plane as indicated in Figure 
D-1A. The one-sigma uncertainty in the S direction divided by 
the hyperbolic excess velociL; produces the one-sigma time of 
flight uncertainty given In Figure D-1B. 
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Figure D-l. Knowledge and Control Uncertainties at Deflection 



The net effect of decreasing the deflection radius on the knowledge 
and control covariances is evident. Decreasing the deflection radius 
from 50- to 30-million kilometers reduces the uncertainties by ap- 
proximately one- third. Decreasing the deflection radius from 30 
to 10 million kilometers produces a decrease in uncertainties of 
approximately one-fourth. 

An intuitive feeling for the exact use of knowledge and control 
covariances may be gained by referring to section B which describes 
the analytical technique used in the Monte Carlo dispersion 
analysis program. 

Execution Errors 


The second source of errors analyzed results from execution errors 
made in implementing the actual maneuvers required. The types 
of maneuvers encountered in the three inodes identified earlier 
are: 

1) Probe deflection maneuver; 

2) Spacecraft deflection maneuver; 

3) Probe release and/or orientation. 

The error models used to analyze each of these maneuvers are de- 
scribed in the following paragraphs. 

The purpose of the probe or spacecraft deflection maneuvers is to 
add a commanded velocity increment Av to che current probe or 
spacecraft state. The actual velocity increment imparted to the 
body will differ from the commanded value by an amount 6Av that 
represents the contribution due to execution errors. The execu- 
tion error model used in this study is defined by four independent 
error sources. 

The first 'rvor source is called the proportionality error and is 
in the direction of the velocity correction, Av , with magnitude 
determined by the proportionality factor, k: 

6Av » k Av 
K 
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A second error source in the direction of Av but independent of 
its magnitude, is the resolution error, s, that corresponds to a 
thrust tailoff error from the thrusters: 


6Av = 
s 


s 

Av 


Av 


The other error sources are pointing direction errors. The first 
of these is a pointing error angle, a, measured in a plane parallel 
to the ecliptic plane (for Jupiter missions, approximately the 
orbital plane) , and along a vector orthogonal to the velocity 
correction vector, Av . If t, J_, k form the unit triad in the 
ecliptic system, then for small angles, a, the velocity error 
caused by the in-plane pointing error is given by 


6Av - c 


A v a Av i- Av i 

L y - x J-J 


- 1/2 


The 


where Av ■ Av 1 + Av 1 + Av k and c * (Av 2 + Av 2 \ 

— x— y z — \x y / 

second pointing error, called the out-of-plane error defines the 

velocity error that is orthogonal to both 6Av and the velocity 


increment vector, Av. Again, for small angles, 3, the velocity 
error resulting from this pointing error, referenced to the 
ecliptic system, is given by 


6Av 


6 



Then the total execution error resulting from a proportionality 
error, k, resolution error, s, and pointing errors, a and 8, is 
given by the sum of these errors 


6Av = 6Av^ + 6Av & + 6Av^ + <5Av^ 

The mathematical model used to describe velocity increment errors 
is the same for either spacecraft or probe. The individual 
magnitudes of the execution error sources, k, s, a, and 8, may be 
varied, however, for individual characteristics of the probe or 
spacecraft. In general, the Pioneer pointing accuracy is 5% of 
the angle rotated off Earth lock, the TOPS error is considered 
to be 1 degree (3a) . 

The probe release and/or orientation error is essentially involved 
with simply aligning the probe axis for its zero degree relative 
angle of attack at entry. This corresponds to the maneuver by 
which the probe is released from the spacecraft in the Mode 3/ 
Deflect Spacecraft sequence or to the probe self-reorientation 
maneuver in the Mode 1/Deflect Probe scheme. The current math- 
ematical model for this type of maneuver is based on a single 
pointing error referenced to the desired direction. 
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Let 11 ^ denote the desired direction for the probe axis. Define 


v 


1 


r x u. 



Then and are unit vectors in the plane normal to the desired 
direction u^. An arbitrary unit vector in that plane may then be 
written. 


x * v^ cos 8 + sin 0 

where 8 is a random variable chosen from a uniform distribution 
over the interval (0, 2 tt). If the orientation pointing error 
is to be of magnitude 6, then the actual pointing direction re- 
sulting from the error is given by 

—A = — D COS ^ * — s * n ^ 


ANALYTICAL TECHNIQUE 


A Monte Carlo technique is used to convert the errors associated 
with the deflection maneuver into dispersions in the critical 
mission parameters. The Monte Carlo procedure will be defined in 
detail for the Mode 1/Deflect Probe deflection scheme. The 
modifications for the Mode 2/Shared Deflection and Mode 3/Deflect 
Spacecraft schemes are then easily explained. 

The Monte Carlo technique use! in this study consists of gener- 
ating a large number of sample probe and spacecraft trajectories 
consistent with assumed statistics of the knowledge and control 
uncertainties and execution errors. These trajectories are then 
propagated to certain mission time points (probe acquisition by 
spacecraft, entry, and selected intermediate points), and the 
critical mission parameters are evaluated. Because each sample 
spacecraft and probe trajectory will differ from the nominal 
(errorless) ones, the critical parameters will also differ. 

The data in dispersions in critical mission parameters are then 
analyzed by empirical formulas, assuming normal distributions 
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to establish mean and three sigma deviations in each critical 
parameter for each selected time point. The computational flow, 
outlined in Figure D-2, may be broken up into three main com- 
ponents that are discussed in detail. 

Preliminary Computations 

The preliminary computations generally comprise the determination 
of the nominal trajectories and preparation for the selection of 
the perturbed trajectories. The nominal deflection maneuver is 
first computed, which for the given hyperbolic excess velocity 
(equivalently the given launch date/arrival date) and desired 
deflection radius satisfies the desired entry conditions and 
communication geometry. The knowledge and control covariances, 
and P £ , are then computed from an orbit determination program 

in which the procedure described in section A is used. Finally, 
the execution error uncertainties, a^, a g , a^, a^, a^, are selected 

on the basis of the spacecraft used, the nominal geometry involved 
(such as the rotation off Earth lock required to implement the 
maneuver), and the magnitude of errors requiring analysis. 

The Sampling Procedure 

The bulk of the dispersion analysis is concerned with the gener- 
ation of the statistically consistent ensemble of spacecraft-probe 
trajectory samples. For each sample, the first problem (listed 
in Figure D-2) is to generate reasonable deviation vectors, 6X^ 

and 6X , from the knowledge and control covariances, P and P . 

C iC c 

This sampling is done as follows. 

Let P represent P^ or P^. Then, since P is positive definite, 
it may be diagonalized. 

P - t t dt 

where D - diag (dj, d|, ... , d|) , T is the orthonormal matrix 
of tie eigenvectors of P, and the superscript T denotes the matrix 
transpose operation. The elements of D are written as squares 
to indicate they are necessarily positive numbers. 



I. PRELIMINARY COMPUTATIONS 

A. Read in nominal spacecraft deflection state X 

nom 

desired entry angle, y, and lead angle, X, 

knowledge and control covariance matrices, P, and P 

k. c * 

execution error model uncertainties, o, , a , a » o*. 

k s or 3 6 

B. Compute nominal deflection maneuver and parameters 


I 


II. SAMPLING PROCEDURE 


Generate large number (N _> 100) of sample off-nominal cases. 
For each case — 

A. Sample knowledge and control covariances to obtain 
the deviation vectors, and 

B. Form the actual and estimated spacecraft states 


act 


est 


X + 6X 
nom c 


X act + 6X k 


Compute the commanded velocity increment, Av c> based on 

the estimated state, % eSt > and the desired y and X. 
Sample the execution errors, generate the resulting 
error in velocity increment, Av. Compute the post- 
deflectior probe state as 


act 


est 


X + 
act 


Av + 6Av 
c 


Y act + 5X k + u 


E. Compute the desired orientation of the probe at entry, u , 


using Y 


Generate a sample orientation error and compute 


the corresponding actual probe orientation, u^. 

Propagate the actual probe and spacecraft to probe entry 
using conic formula. Store the dispersions in critical 
mission parameters at each timepoint. 


III. STATISTICAL ANALYSIS OF DISPERSION DATA 

Using the dispersion data generated for each sample case, 
compute the mean and standard deviaiton for each critical 
mission parameter for each selected timepoint. 


Figure D-2 Computational Flow Chart of Monte Carlo Analysis Program 





T defines a transformation from the original Cartesian coordinate 
system to a new system in which the covariance matrix is uncor- 
related, thereby allowing the individual components to be sampled 
independently. A vector random variable from a distribution of 
mean zero and covariance, D, is given by 

Z - (ei di, 62 d2, ••• , e$ dg) T 


where each e^ is a scalar random variable sampled from a normal 
distribution of mean zero and standard deviation unity. The e^ 
are computed from the formula 

e i 

where and are random numbers generated from a uniform dis- 
tribution over the interval (0, 1). The correlated deviation 
vector in the original Cartesian coordinate system may now be 
computed using the transformation matrix, T, as 


^-2 in cos g^ 


6X - T6Z 


Therefore, for each sample case, the knowledge and control 

covariances may be sampled as above to produce knowledge and 

control state deviation vectors, 6X, and 6X . The "actual" 

k c 

spacecraft state for the sample is given by 


X 


act 


X 

nom 


+ <$X . 
c 


The "estimated" state of the spacecraft at deflection is given by 


X t - X + 
est act K 


Furthermore , 
the ensemble 


fct the end of a large number of samples 
of deviations, 6X c and 6X^, should obey 


(N > 100) , 
the empirical 


formulas 


E |<5X c 6X c T j 

E <) 


-+ P as N -*• » 
c 

P. as N + ®. 
k 


These approximations are checked at the end of the sampling pro- 
cess to ensure that a statistically consistent set of data has 
been generated. 
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Having obtained the estimated state of the spacecraft for a sample, 
the next task is to determine the commanded velocity increment, 

Av ^ t to be added to the probe at deflection. The vector, Av c> is 

computed as a function of the estimated spacecraft state at de- 
flection, the desired probe flight path angle at entry, y # and 

the desired lead angle, A c . The lead angle is the angle between 

the radius vector to the probe at entry and the radius vector to 
the spacecraft at the time of entry. The lead angle is negative 
when the probe leads the spacecraft. The deflection velocity is 
determined by iterating on the value of A0, the true anomaly 

increment of the probe in going from deflection to entry. The 
true anomaly of the probe at entry, f , can be expressed as: 

rh 


PE 


V 


r ej* Ae * r a 


is the radius of the bus at deflection and R. is the 

A 

radius of the atmosphere, f R , A9, R are used to compute 

F E , EJ A 

the time of flight, T, from deflection to probe entry using the 
universal form of Kepler’s equation as presented in Reference 1. 
The spacecraft state is then propagated forward for time T, again 
using the universal form of Kepler’s equation, and the angle, X, 
is computed from the state of the spacecraft and the state of the 
probe at time T. The angle, X, is compared with the desired X^ 

and the A0 which causes X to be within 0.01 degrees of X^, is 

found. Since the orbital elements of the probe are known, the 
required deflection velocity vector may be calculated. The com- 
manded deflection velocity is in the spacecraft plane of motion. 

The deviations in the implement'd deflection velocity from the 
commanded are generated using the model described in subsection 
A2. This model assumes the implementation error is given by 
three sources: proportionality error, K; resolution error, S; and 
pointing angle errors, a, 3. These errors are assumed uncorre- 
lated and normally distributed. The standard deviation of each 
error o , o f o , a is input to the program. The error used 

K b Ot p 

in each Monte Carlo sample is generated by 
K, S, a, 6 * ei e 2 o g , e 3 0 a , 

where each is a scalar random variable sampled from a normal 

distribution of mean zero and standard deviation unity. The 
execution error, 6Av, is then the sum of these errors. 


where 


"ej 
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The actual state of the probe at deflection is then the actual 
position vector of the probe and the sum of velocity of the 
spacecraft and the actual deflection, Av, given to the probe. 


act 


X + 
act 


j^Av c + 5Av 


Now that the actual states of the probe and spacecraft » are known, 
the time histories of these two point masses may be computed 
using conic models. 

The attitude of the probe is a crucial parameter, however, as this 
determines such mission critical parameters as angle of attack 
at entry and probe aspect angle. Using the nominal trajectory, 
the desired probe attitude is computed, u^. A pointing error 

caused by imperfections in the attitude control system is computed 
in a similar statistical fashion as described above. Once the 
"actual" probe axis is computed, it is assumed to hold that 
orientation throughout the mission. No errors are added to the 
spacecraft axis since it is assumed that the spacecraft can hold 
the Earth lock with essentially no error. 

For each sample probe and space trajectory and probe attitude, 
the critical entry and communication parameters may be computed 
as functions of time. The resulting collection of data must 
then be analyzed by the techniques described in the next sub- 
section. 


Statistical Analysis of Dispersion Data 


The empirical computation of the standard deviations of scalar 
parameters such as entry angle, lead angle, or lead time may be 
computed in the following way. Let a. be the value of any such 

th * 

parameter for the i case. Suppose that there are N samples 
to be analyzed. Then the mean and standard deviation of the 
distribution of a are given by 


•-* E ? i 

i - 1 


N 


y* V 

N 


i •* 1 


_ 2 
- a 

N 
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These formulas must be extended for vector quantities. The probe 
entry site dispersions are given in terms of the two vectors of 
latitude and longitude (LAT, LON). The spacecraft- probe look 
direction is conveniently described in terms of the two vectors 
of cone angle and clock angle (CA, CLA) referenced to Earth and 
Canopus. Equivalently, these dispersions may be defined in terms 
of cone angle and cross-cone angle (CA, CCA). 


Let Z represent the vector describing the actual values achieved 

for a vector quantity on the i sample- Then, the vector of mean 
values and the covariance matrix describing uncertainties and cor- 
relations of the vector are given by the formulas 


2 . i E 

N 


i - 1 


P - E(ZZ T ) - E(Z)E(Z) T 

z r n i 



These formulas are used to compute the covariance matrices of the 
critical mission vector quantities as well as to reconstruct the 
original knowledge and control covariances from the deviation vec- 
tors generated. 

For the entry stie dispersions (LON, LAT) and the spacecraft-probe 
look directions (CA, CLA) or (CA, CCA), the 2*2 covariances are 
further analyzed. Let any such covariance matrix be denoted P. 
Then since P is positive definite it may be diagonalized to pro- 
duce 



Where and X„ are the (positive) eigenvalues of the covariance 

A D 

matrix P, and T is the orthonormal matrix of the eigenvectors of 
P. Let 
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Then t A Is che unit eigenvector associated with the eigenvalue, X A * 

Then the angle between the vector, t A « and the u x unit vector 

(longitude direction for (LON, LAT) , cone angle direction for (CA, 
CLA) or (CA, CCA) is defined by 

6 ■ cos” 1 

The uncertainty ellipse may now be easily constructed: X A repre- 
sents one semiaxis i X fi the other. 


) 0 <_ 6 _< it 



Pigux-e D-3 Tuo-Dimensional Uncertainty Ellipse 


Procedures for Alternative Deflection Modes 

The procedure for Mode 1/Deflect Probe has been described in detail. 
A comparison of the procedures for the Mode 2/Shared Deflection and 
Mode 3/Deflect Spacecraft is provided in Figure IV- 3 of Volume II 
of this report. 

For each sample in each mode, the initial task is to determine the 
control and knowledge deviations, <5X c end 6X^, by sampling the con- 
trol and knowledge covariances. The actual state of the spacecrart 

is then defined by X „ - X + 6X ; the estimated state of the 

act nom c 

spacecraft is X - X _ + 6X. , where X is the nominal space- 

r est act k nom 

craft deflection state. 
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In Mode 1/Deflect Probe the estimated state, X , is used to de- 
termine the commanded probe deflection, Av. This velocity increment 
is then degraded by an execution error, <SAv, determined by sampling 
from possible execution errors. The probe axis orientation through- 
out the mission is computed as the nominal orientation corrupted 
by an orientation pointing error, 6. The deflection states for the 
probe and spacecraft are then propagated to a series of time points 
at which the sample deviations are recorded. These dispersions are 
then analyzed to yield mean and standard deviation values. 

In the study of the Mode 1/Deflect Probe deflection scheme, it was 
discovered that the knowledge and control uncertainties did not 
significantly affect the critical mission parameters. Therefore, 
to simplify the Mode 2/Shared Deflection analysis, it was assumed 
that the knowledge and control uncertainties were zero so that the 
spacecraft correction, Av, could be precomputed and simply read in. 
At probe deflection, the probe is commanded to be aligned in the 
direction of both the Av and zero angle of attack orientation. 
Because of the in-plane and out-of-plane pointing errors, the cor- 
rect direction is not achieved. The incorrect orientation is 
then used for both the Av addition and the probe longitudinal axis. 
Following the probe deflection, the spacecraft correction is im- 
plemented with the execution error model described in subsection 
A2. The resulting sample deviations are then collected and anal- 
yzed to determine the mean and standard deviations. 

For Mode 3/Deflect Spacecraft deflection analysis, the knowledge 
and control uncertainties are again considered. The probe is ori- 
ented at release using the input orientation pointing error, 6, 
and keeps that attitude throughout the mission. The estimated 
state of the spacecraft, X _ gt » * s used to compute the deflection, 

Av. Spacecraft execution errors are then sampled to determine 
the error, 6Av. The erroneous velocity increment is then added 
to the spacecraft and the spacecraft is then propagated to the 
selected time points. The critical mission parameter deviations 
are then recorded and analyzed to determine the important 
statistical data. 
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TURBOPAUSE DESIGN MISSION 1A, RUN 5 (ACTUAL) 

PHASE 6 ENTRY ANGLE 23.2U5 DfG INERTIAL AT 72550. ? KM RADIUS 
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SCALE HEIGHTS (NEUTRALS) HI -54.974 H2 -27.487 HE -13.787 
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INSTRUMENTS LANS PRUDE MASS SPECT IRPA NRP A OPT SPECTk 
KM/MEASURFMERT 11.445 9.159 11.449 11.449 68.694 
NO MEAS TOTAL 130.254 2.163 112.647 45.312 125.73U 
NO MEAS PHASE 1.746 2.163 1.746 1.746 .251 



TURB0PAU3E OESI&M MISSION 1A, RUN * <ACTUAtl 
PHASE 7 ENTRY ANGLE 2*. *>05 OEu INERTIAL AT 72GG0.0 KH RADIUS 
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TURBOPAUSE OESIGN MISSION 7, RUN 6 CACTUAL) 

PHASE 2 ENTRY ANGLE 33.300 OEG INERTIAL AT 7255C.0 K.M RADIUS 
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ALTITUDE 1800.0 KM T IMP .3731 HOURS TI*E .C135 HOURS TIME .0017 HOURS V I NE RT I Al_ 59.8401 KM/SEC 

RADIUS 7335C.0 KM FROM 22.388 MIN TO .807 MIN IN .100 MIN VREL AT I VE 49.5807 KM/SEC 

RX<1> 1343.275 SEC TA*G 48.44? SEC PHASE 6.003 SEC VPADIAL 33.2856 <*/SEC 
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ALTITUDE 24 0.0 KM TIME .3«62 HOURS T I*E .0003 HOURS TIME . 00 05 HOURS VINERTHl 60.4731 KM/SEC 

RADIUS 71790.0 KM FROM 23.175 MIN TO .020 MIN N .030 MIN ^RELATIVE 50.2623 KM/SEC 

RX<11 1390.496 S£C TA»G 1.221 SEC PHASE 1.030 SEC VRAPIAL 32.7763 KM/SEC 
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TUR90PAUSE DESIGN MISSION 7, PUN 6 (ACTUAL) 

PHASE 6 ENTRY ANGLE 33.300 OEG INERTIAL AT 7?5S0.0 KN RADIUS 
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TUR90PAUSE Of SIGN MISSION 7 # RUN 6 CACTUAO 
PHASE 8 ENTRY ANGLE 33.300 OEG INERTIAL AT 72550.0 KM RADIUS 
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APPENDIX F 


GSFC ATMOSPHERIC MODEL 
(ATTACHMENT II TO STATEMENT OF WORK) 


January 4, 1971 


Altitude ab 



Temperature, °K 


Fig, F-l Temperature of Jovian Atmosphereo 



Altitude above Turbopause 



Fig. F-2 Molecular Weight of Jovian Atmospheres 
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Fig. F-3 Number Densities of Jovian Atmosphen 




Fig, F-4 Number Densities of Jovian Atmosphere 
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%■ F ~5 Number Densities of Jovian Atmosphere 
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Log of Number Density, 1/cm 
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Fig. F-7 Number Densities of Jovian Atmosphere 



Fig . F-8 Niwiber Densities of Jovian Atmosphere 






Fig. F-3 Density and Pressure of Jovian Atmosphere 
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Fig. F-10 Density and Pressure of Jovian Atmosphere 
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Fig . F-ll Density and Pressure of Jovian Atmosphere 
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PROBE K-BAND POWER SOURCES 


R. Compton 


January 6, 1972 


Preliminary analysis of the atmosphere entry communication black- 
out problem for the Jupiter Turbopause Probe indicated that fre- 
quencies in the K-bana region may be required to provide a data 
link below the turbopause withcut excessive (> 3dB) attenuation. 
Lower frequencies will have greater attenuation since plasma at- 
tenuation is inversely proportional to the frequency of operation. 
Therefore, a data transmission system operating at K-band would 
be necessary if deeper probe penetration were required on some 
future missions. 

Vendor Survey 

A vendor survey was performed to determine the projected 1975 
state-of-the-art for RF power sources in the centimeter region. 
This survey was necessary to determine the usable upper RF power 
limit in designing the missions. Both so.-id state and tube 
sources were investigated. General industry trends were sur- 
veyed along with the feasibility of achieving by 1975 a trans- 
mitter with performance characteristics as listed in Table G-l. 
K-band was chosen Decause it was the upper frequency considered. 
Eight suppliers were contacted and six responded to the request. 
Summaries of their responses follow. 

Sperry Electronic T u be Division, Gainesville, Florida - Sperry 
is presently engaged in development of a traveling wave tube 
(TWT) to achieve an objective of 50 watts cw from 18.0 to 26.0 
GHz. It will be designed for either conduction or air cooling. 

The gain at rated power will be 30 dB minimum. Sperry antici- 
pates that a demonstration model will be complete by April 1972. 
The TWT will be periodic permanent magnet focused. Dimensions 
and weight have not been established but would probably be within 
12 x 2 x 2 in. and 5 lb, respectively. The TWT would be designed 
to serve both ground based and airborne applications. The tube 
will have a helix interaction circuit to provide the required 
bandwidth, and is expected to exhibit a minimum of 20% overall 
efficiency at bandedge power levels rising to over 25% at center 
frequencies . 

The phase and intermodulation characteristics of this tube should 
be similar to those exhibited by Sperry’s other high power TWTs 
at lower frequency ranges. Minimum power consumption and weight 
are requirements. In narrower frequency ranges of 1 GHz and less, 
several techniques such as helix phase tapering and multiple col- 
lector depression can raise TWT efficiency to near 40%. However, 
for greater bandwidths, it would be very speculative to assume 
more than 30% TW^ efficiency by 1978. Power supply efficiency of 
85 to 92% would put overall transmitter efficiency somewhere in 
the 25 to 28% region. 
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Table G-l. Transmitter Performance Characteristics 


PARAMETERS 

Frequency 

RF power output 

Modulation 

Efficiency 

Primary power 

Size and weight 

Life 


K-band (18-25 GHz) 

25 watts cw 
Phase 

Maximum feasible 
18-26 V dc 
Minimum feasible 
18 months nonoperating during 
outer space cruise; 48 hours 
operating over 15-day period 


ENVIRONMENT 

Ambient pressure 
Shock 

Sine vibration 


Random vibration 

Thermal vacuum 

Decompression 
Temperature change rate 
Acoustic 


10"“ 14 torr 

2Q0g, 1 ms 

5-35 Hz, 0.75 g 

e rms 

35-250 Hz, 6 g 

e rms 

250-600 Hz, 4.5 g 

°rms 

600-2000 Hz, 9.0 g 

& rms 

(3 axes, 1 octave per minute) 

0.3 g 2 /Hz, overall 16 g 

60 sec each axis, 3 axes 

-20°C, 10“ 6 torr, 2 days; 

1 day at +55°C 

124 to7;r/sec 

30°C/hr 

150 dB for 2 minutes 
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Sperry f s development plan in solid state microwave devices includes 
wideband amplifiers to serve many of the applications that now 
require traveling wave tubes. Power levels of 20 watts at K-band 
in six to eight years would not seem to be impossible, but they 
cannot comment with any assurance on efficiency and the other 
characteristics that would be vital to our application. 

General Electric Company, Owensboro, Kentucky - GE has made a re- 
cent study of the available and theoretically available perfor- 
mance from solid state sources, and 20 watts appears to be beyond 
the present state of the art. The 1975 predicted singlechip 
performance was in the 5-watt area, which is still short of the 
requirements. The industry proposes 4 watts at 12 GHz which 
should be available soon. This is a Fairchild development. It 
uses parallel and/or series operation to obtain the power output. 
This suggests that 20 watts may be available in the future from 
a composite source using hybrid power "trees" to couple several 
diodes Lo one output port. 

EMR Telemetry, Sarasota, Florida - A search of microwave litera- 
ture indicates that technology is available to produce 20 watts 
cw in the band of 18 to 25 GHz at the present time. This can 
be accomplished by use of a power oscillator operating at the 
desired output frequency, or by generating a low power signal 
which is then used to drive an amplifier that has the gain and 
desired power output capability at the output frequency. 

In either of these cases, electron tubes would be required. The 
power oscillator might be either a backward wave oscillator or 
a multicavity klystron, and the power amplifier could be either 
a traveling wave tube or a multicavity klystron. 

At this time there appears to be no practicable way that 20 watts 
of cw power can be produced over this frequency range by solid 
state methods. Present varactor diodes that operate at these 
frequencies have a very low reverse voltage capability; therefore 
they have a very low power capability. Solid state diodes that 
can generate power directly at the frequency are use severely 
power limited. 

A breakthrough in the solid state field may make powers in the 
order of 20 watts at K-band available by 1978, but at the moment 
there is no way of predicting the characteristics of such a 
system. For this reason, our predictions will be made on the 
basis of further development of electron tubes for the production 
of high power at K-band. 
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The projected transmitter would use a diode oscillator to produce 
a low level signal at the output frequency. This low level signal 
would be stabilized by sending part of the signal to a discrimin- 
ator, also operating at the output frequency, where any frequency 
deviation would produce an error signal. Tnis error signal would 
be used to keep the output frequency within the desired tolerance. 

The low level signal would be amplified to the desired level of 20 
watts in a TWT amplifier. Phase modulation could be accomplished 
by the variation of one of the potentials applied to the TWT. It 
is estimated that the overall efficiency of such a transmitter 
would be in the order of 8 to 10%, the primary power required 
would be about 250 watts, and the size. 0.15 ft 3 . The weight would 
be 10 lb maximum. 

Except for the 18 months in space, the environmental requirements 
are not stringent. The best means of ensuring operation after 18 
months in space might require a slow leak in the transmitter hous- 
ing to permit the evacuation of the housing during the 18 months 
without the danger of refilling with Jupiter’s atmosphere during 
entry . 

EMR is developing a solid state S-band transmitter in which the 
power is developed by an oscillator operating at the output fre- 
quency. This oscillator is stabilized by sending p_rt of the os- 
cillator signal to a discriminator which operates at a harmonic 
of the output frequency. The error signal developed in the dis- 
criminator is used in a closed loop system to reduce the undesired 
frequency deviation of the S-band oscillator. The experience 
gained in this development program should be helpful in develop- 
ing a closed loop control system for a K-band oscillator. 

Hughes Electron Dynamics Division, Torrance, California - Hughes 
Electron Dynamics Division presently manufactures traveling wave 
tubes that supply 20 watts at 14 GHz. This tube is designated 
the Hughes Model 851H and is used in the Skylab Program. Hughes 
also has traveling wave tubes that deliver 10 watts at 18 GHz. 

Tiiey see no problem in extending performance to 20 watts at 18 to 
25 GHz (Figure G-l) with high reliability. 

In the past several years, efficiency achievable in a TWT has 
increased in the range of 35% to a meacured 60% at narrow band. 

High efficiency at K-band is still difficult because of high cir- 
cuit losses and the small size of the helix. Hughe- expects to 
achieve 50% efficiency in 1972 over a band centered at 6 GHz and, 
based on this, can confidently project an efficiency of 40% at K- 
band by 1978 to 1980. The size will be determined by the required 
gain; however, for a 40-dB gain tube, weight would be approximately 
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1.25 lb. Power supply efficiency would be in the order of 85 to 
90%. The complete package, tube and power supply would weigh ap- 
proximately 10 lb and would measure 3.5 * 5.0 * 10 in. The length 
would be dependent upon gain of the tube. Regarding the environ- 
mental requirements, present hardware will meet these specifications. 

Conic Corporation, San Diego, California - By using a Gunn source 
and hybrid combining techniques of Impatt amplifiers, it is Conic* s 
opinion that 20 watts will be achievable using solid state tech- 
niques in the 1978 to 80 time frame. It is also quite possible 
that there will be Trapatt diodes available for high power gener- 
ation at K-band. At present, the GaAs Impatt diode shows the 
greatest promise for high cw output. The achievable efficiency 
will be in the 15 to 20% region. Primary power requirements will 
still be in the 60 to 90 V dc region. Strip line circuit tech- 
niques will reduce overall weight to 7 to 10 lb and occupy 100 to 
150 in^. The environmental requirements are no more severe than 
those imposed on present day transmitters. 

Watkins-Johnson Company, Palo Alto, California - Watkins-Johnson 
stated that the requirements listed in Table G-l could be met by 
a TWT amplifier using present design techniques. They presently 
have a line of TWTs (WJ-231 series) in development that have a 
saturated power output of 35 watts cw at 8 GHz(X-band). The over- 
all efficiency is 28% and the weight is 2.5 lb. The noise figure 
is 30 dB maximum and has a VSWR of 1.4 to one. 

For K-band the projected design details include an overall effi- 
ciency of 30%, primary power input for 20 watts output is 68 watts, 
overall weight (amplifier + power supply) is 7.5 lb, and will meet 
all the environmental requirements. 

Literature Survey 

A literature survey was also conducted to obtain a more objective 
opinion on the projected state of the art for K-band RF power 
sources in the 1975 to 1978 time frame. Design approaches are 
divided into two categories; microwave tubt ? and solid state de- 
vices. Microwave tubes include only two devices; the solid state 
category is more general and includes a broad range of devic s. 
Performance characteristics of the principal RF sources are shown 
in Figure G-2 and will be discussed in detail. 1 




a i co *h 
^ C V) 
O flj CO 

u 

eui 



nj 

T3 CO 
T3 0) 
•H 

O H 



Projected 1975 RF Power Source Capability 













Microwave Tubes - The reduced performance of standard grid-con- 
trolled vacuum tubes in the microwave region has fostered devel- 
opment of other types of tubes for use as amplifiers. Such tubes 
function on the basis of the modulation of the velocity of an 
electron stream rather than on its density. Microwave tubes are 
divided into linear beam and crossed-f ield devices. The princi- 
pal types of linear beam tubes are the klystron, the TWT, and 
the backward-wave oscillator (BWO) . Crossed-f ield tubes include 
the magnetron oscillator and carcinotron oscillator. The magne- 
tron requires a strong magnetic field to develop the required 
power which would interfere with the science instruments on the 
probe and therefore was not considered further. Carcinotrons are 
basically a laboratory device that have an inherently high noise 
level with spurious power output, often "Ot more than 10 to 15 dB 
below the main signal output.^ Therefore, carcinotrons were not 
considered further for probe transmitter applications. The BWO 
is generally a low-power device, with 50 mw being typical; it is 
not suitable for high-power applications. 

Of the several types of microwave tubes, only the klystron and 
TWT are possible candidates for the transmitter. Multicavity 
klystron amplifiers find considerable application for high-power 
cw power systems. Average powers of 1 watt or greater are pre- 
sently available at X-band with an efficiency of 40%. Klystrons 
have the disadvantage of restricted bandwidth, high beam voltage 
(enhances RF breakdown during atmospheric ^ntry) , and frequency 
instability. The 20-watt klystrons must also be liquid cooled 
for best results and are electrostatically focused. 

A TWT differs from a klystron in that the RF field is not confined 
to a cavity but is distributed along a wave-propagating structure. 
An electron gun produces a longitudinal beam of electrons that 
interacts with the field of a wave traveling along the structure. 
Spaced closely around the beam is a helix that propagates slow 
wave (Phase velocity is appreciably less than the velocity of 
light )• An external RF signal is coupled into the tube where its 
longitudinal component interacts with the electron beam. Mutual 
interaction exists between the electron stream and the wave with 
the net result that dc energy is released by the electron stream 
to the circuit as RF energy, and the wave is thus amplified. TWTs 
are designed to operate in the power saturation region since effi- 
ciency is highest at that point. Raising the input power above 
the saturation point will lower the RF output power and effi- 
ciency. 4 , 5 
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TWTs have a broader bandwidth than klystrons, typically 1 GHz, 
and have a phase linearity of +3 deg over any 10-MHz bandwidth 
within the frequency range of the tube. A TWT has a longer life 
span than a klystron. Presently, TWT efficiency is slightly 
lower than for a klystron, but TWT suppliers are actively work- 
ing on raising tube efficiency to levels of 40% or greater. Much 
of this impetus stems from the growing need for high power tubes 
for space applications. 

Much of the power wasted in a TWT is in the form of heat gener- 
ated by the impact of beam electrons with the collector structure. 
Velocity tapering is used at the collector to reduce electron 
impact velocities through multistage collector depression (MSCD) . 
One of the advantages of collector depression is that it lowers 
the demands on high voltage power supplies. Heat dissipa .ion in 
space is generally accomplished by conduction and radiation. The 
use of heat pipes for thermal control in a space vehicle will be 
tried for the first time in a tube being developed by Hughes. 8 

Hardware development activity is increasing to develop a 12-GHz 
tube producing about 20 watts for the newly authorized K-band 
channel for Comsat and Intelsat-5. It is interesting to note 
that all tube suppliers are proposing TWTs for this application 
in lieu of klystrons. 7 This coincides with the recommendations 
developed from the vendor survey conducted for the probe study 
considering TWTs. Klystrons or solid state devices were not 
recommended by any supplier. 

Solid State Devices - Two trends are putting solid state devices 
into microwave circuits in increasing numbers: (1) the increased 

use of integrated circuits leading to smaller, lighter s and more 
rugged equipment; (2) the increased capability of solid state 
devices to perform at these frequencies. 8 As mentioned in the 
vendor summaries, the state of the art for solid state devices 
for 1975 is not predicted to have high power RF amplifiers 
available at K-band as reliable, space qualified production items. 
Industry emphasis is on the more popular communication bands (X 
and S) which are lower in frequency. 8 Several devices that show 
promise of operation in the microwave region exist. These devices 
will be summarized briefly. 

Diodes that operate under the tra' r rred-electr on bulk effect 
using GaAs operate up to 35 GHz, v presently only have power 
outputs of 0.5 watt or less. Gunn-effect devices have their best 
performance at higher frequencies than do transistors. The 
extrapolated 1975 cw Gunn-effect capability is presented in Fig- 
ure G-2. 10 
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Semiconductors operating at microwave frequencies are transit- 
time limited. With single-chip power silicon transistors, 5 
watts at 2 GHz have been obtained. Germanium planar-signal 
transistors are also being developed. Conversion efficiency 
is a critical factor. Projected 1975 data on power versus fre- 
quency is plotted in Figure G-2. X-band is the projected 
upper limit with an efficiency of 40% and producing 1 or 2 watts. 

Limited space charge accumulation (LSA) devices of bulk GaAs 
oscillate because they exhibit a negative-resistance character- 
istic and are not transit-time limited. They reach a peak in 
both output and efficiency above 8 GHz with an upper limit at 
500 GHz. They are presently a low-power device and have pro- 
duced a peak power of 50 watts at X-band. The average cw power 
is 5 watts. Figure G-2 depicts the oerformance of LSA devices 
for 1975. n , 12 

Another device that exhibits a negative resistance and avalanche 
breakdown is the impact avalanche transit time (Irapatt) diode. 

The fundamental limit for these devices is about 20 watts at K- 
band, which has not been reached to date. Efficiencies as high 
as 25% appear possible with a 1975 projected efficiency of 20%, 
as seen in Figure G-2. 13 , 14 

Since the anomalous mode or trapped plasma avalanche triggered 
transit (Trapatt) mode of avalanche diode operation was dis- 
covered only about three years ago, sufficient trends have not 
been established to accurately extrapolate performance to 1975. 
Progress is being made in extending the frequency of operation of 
such devices. 2.4 watts cw were generated at S-band in 1970 with 
an efficiency of 20% using silicon. Theoretical efficiencies are 
in the region of 50 to 60%. RCA has predicted a capability at 
10 GHz of 100 watts cw with 50% efficiency. This device may offe*" 
the best means of competing with tubes on a realistic cost per 
watt basis during the next decade. i2 , 15 

A combination of two solid state devices also offers possibilities. 
A Gunn diode master oscillator, frequency locked and modulated by 
a varactor diode, followed by an Impatt diode to amplify the signal 
could be used. Approximately 5 watts is projected for the 1975 
time period at about 15 to 25% efficiency per Impatt diode. Sev- 
eral devices can be cascaded to provide a higher total RF power, 
in the normal failure mode, an Impatt diode produces a momentary 
short circuit (a few milliseconds) followed by an open circuit. 
Therefore, loss of any one diode will result only in the loss of 
power contributed by that circuit element. 
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Conclusion 


As can be seen from vendor responses and the literature survey, 
the best candidate for a K-band RF power amplifier in the 20- to 
25-watt region is a traveling wave tube amplifier (Figure G-2). 

As seen from the illustration, 25 watts at K-band does not exceed 
the tube state of the art, and the 1975 projected solid state 
capability may also meet the power requirements. 
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The inviscid body flow field is defined as that region bounded 
by the bow shockwave and the body boundary layer. As noted in 
Figure H-l, the flow field is composed of a subsonic region and 
a supersonic region. In the nose cap (subsonic) flow, the govern- 
ing partial differential equations are elliptic, and in the after-’ 
body (supersonic) region they are hyperbolic. Therefore, dif- 
ferent computational techniques are used to describe the two flow 
fields with suitable matching procedures applied at the sonic 
line. The details of this computational method are presented in 
this appendix, following a description of the nonequilibrium 
thermochemistry model. 


THE THERMOCHEMICAL MODEL 


The gas in the inviscid shock layer consists of monatomic and 
' diatomic species, including ions and electrons. For this analy- 
sis. the species are H, He, e, H 2 , H + , and He + ; however, the 
method to be described could be applied to any such gas mixture. 
Each component of the gas is assumed to be thermally perfect. 
Thus, the equation of state for the gas mixture may be written 


N 

s 

P = p RT £ Yj [H-1J 


where R = 8.3134 x 10 7 ergs/mole-°K is the universal gas constant, 
Yj - the species molal concentration in mole/gm, p is the gas 


density, and P and T are the gas pressure and temperature, re- 
spectively. The caloric equation is written in the form 


h. - h? + e. + I. [H-2] 

J J J J 


where e. is the vibrational energy, h is the heat of formation, 

3 3 

and Ij represents the energies of translation, rotation, and elec- 
tronic exitation. Therefore the 1^ term expands to 
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where denotes the number of atoms per molecule of the j ^ 

th 

species, e is the energy of the £ electronic level, g is 

the degeneracy of this level, and m is the number of electronic 
levels being considered. 


The chemical potential 
equation 
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of the j species is given by the 
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The constant, b ^ , includes the translational and rotational con- 
tributions, in addition to the logarithm of the ground-state de- 
generacy, and may be expressed as. 


bj = £n 


2.56033 x 10 2 m 3 / 2 g _ /% 
j B jl/ r.j 


"j ' 1 


, th 


where m. is the molecular weight of the j species, and 6 is 

J th ^ 

the characteristic rotational temperature of the j species. 
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For species in vibrational equilibrium, the simple harmonic oscil- 
lator approximation , 


£ . 
J 




[H-6 ] 


th 

is used for the calculation of vibrational energy of the j 
species. Note that 0 is the characteristic vibrational tempera- 

th Vj 

ture of the j species. The rate of change of concentration of 

the j th chemical species along a given streamline, s, in the flow 
field is governed by two equations. 
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X jk 3s 


0 , 


j = N + 1, N +2, .... N 

c c s [H-7J 


k = 1, 2, . . . , N £ [H-8] 


where N^ represents the number 

the number of atoms of the k th 

N is the number of chemical reactions. Conservation of electric 

R th 

charge is assured by assigning a., = -1 if the j species is a 

^ e 

positive ion. Note that the convention for ordering chemical 
species in the model lists the monatomic elements first, followed 
by the molecular species and ions. Note also that in this order- 
ing, the electron is considered an elemental species. 

1. Chemical Reactions in the Flow Field 


of elemental species and a_.^_ is 
th 

element in the j species, and 


The set of chemical reactions may be represented by the expression 
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• V V. .* M. , 

3 J u t-* ij J 


b, j=l 


i = 1, 2, ...» N r [H— 9 ] 
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where Mj represents the j 


th 


chemical species in the reaction, v 


ij 


and v^* are respectively the stoichiometric coefficients of the 
j ^-reactant and the j ^-product in the i^ 1 reaction, and k and 

th 1 

k^ are the forward and backward reaction rates for the j re- 


Note that N D is the number of chemical reactions in the 
K. 


action 

system, and j is the number of chemical species being considered 


The forward and backward rate constants are related by 



[H-10] 


where AF^ is the chanr Ln free energy in the i^ reaction. The 

subscript "eq" on the reaction rates indicates that the rates are 
evaluated for equilibrium vibration. Because the finite rate of 
relaxation delays the population of the upper vibrational levels , 
a greater effective energy per collision is needed for dissocia- 
tion. Hence, the equilibrium rate constants must be modified by 
a coupling factor. Therefore, 


k 


f 


i 




[H-ll ] 


where the exponent, A , may be either zero or one to indicate 

whether or not the given reaction is affected by the vibration- 
dissociation coupling factor, V ^ . By using Equation [H-10] the 

forward or backward reaction rate constant may be specified (in- 
put) and the other reaction rate computed from it. For example, 
the forward reaction rate constant can be expressed as 


B . 


V 


exp 


- [ci/RT D1 J 


eq 


[ H— 12 ] 
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where is called the frequency factor, is the activation 
energy for the reaction, and the temperature exponents, B. and 
, are usually obtained from correlation of experimental re- 
action kinetic data. With these results and Equation [H-7], the 
rate of production of a given species along a given flow field 
streamline, s, may be more explicitly written as 



In Equation [H-13], the parameters, W^, D^., and Z ^ may be either 

zero or one to denote the type of reaction. For example, in the 
case of specific collision partners, = 1, while D = Z « 0. 
Such a reaction would be of the type 


[H-13] 


H 2 + H £ 2H + H. 


If all the collision partners are lumped together, such as in 
H 2 + M t 2H + M, 


then Dj = 1, while W = = 0 


The third term is used when some of the collision partners are 
lumped, while others are specified. 
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The parameter appearing in Equation [H-13] is defined as 
N 

— nw ’ 11 " 11 

r 

j 

J-l 

where is the equilibrium constant for the i ul4 reaction. The 

significance of is that it represents the degree of nonequilib- 
rium, 

= l - averse rate 
x i forward rate 

and at equilibrium, is zero. Immediately behind the shock, 

X i = 1; in a region of rapid recombination, ma y assume large 
negative values. 


ij 




n-1 



B. THE FLOW FIELD EQUATIONS 


As noted earlier, the inviscid body flow field is defined as that 
portion of the flow bounded by the bow shockwave and the body 
boundary layer. The flow forward of the sonic line (Fig. H-l) is 
subsonic; that aft of the sonic line is supersonic. As a result, 
the partial differential equations governing the flow field change 
from elliptic to hyperbolic as the flow crosses the sonic line. 
Since this is the case, the analysis is best performed by using 
different integration techniques in the two regimes. In the sub- 
sonic regime the equations are numerically integrated using 
Runge-Kutta techniques; for the supersonic flow a modification 
of the method of characteristics is applied. These two flow re- 
gimes are discussed in this section. 
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Bow Shockwave 


Sonic Line 


Boundary Layer 



Figure ti-1 Inviscid Body Flow Field 

The Inviscid Subsonic Body Flow Field 

The inverse method, which specifies the shock shape and calculates 
the body shape, was selected for this analysis primarily because, 
in contrast to the integral method of Belotservkovskii, the inverse 
method yields detailed results at field points throughout the shock 
layer. The differential equations governing the flow field are 
integrated on a coot Lnate system in which one set of rays is 
parallel to the shocx and the other set of rays is normal to the 
shock. The thermochemical flow field properties rre then calcu- 
lated at each grid-point of this curvilinear, orthogonal coordi- 
nate system. Figure H-2 illustrates the coordinate system, and 
Figure H-3 shows the variation of temperature and density at 
selected grid points in the subsonic flow field. 

The vector forms of the subsonic flow field conservation equations 
follow. 

Conservation of Mass 

V • (p'u) = 0 £ H— 15 3 

Conservation of Momentum 

(u > V)u = ^ VP [h-16 ] 

Conservation of Energy 

u • VH = 0 [H-17 ] 
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Figure H-2 Coordinate System for Integration of the Subsonic 
Flow Field Equations 



Figure H-S Shock Layer Field Point Coordinate System with 
Calculated Thermo-Properties Variations 


- 95.5°K 
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Conservation of the j 


Chemical Species 


th 


" • w ■ p r Q u ( t> o> hjj- c j) 

i»l r 

Conservation of Chemical Element 
N 

s 

£ « jk (u ■ V Tj ) - 0, 

j-1 

Conservation of Vibrational Energy 
/ '■ 

u ‘ ' c 3 * ”3 r> c 3’ L Q 


ij 


i-1 


k - 1, 2, 


The two equations of state supplementing this system are, 
T 

and 


= A(P) (MW) 

P 


[H-18] 


, N c [ H— 19 ] 


[H-20] 


[H-21] 


N 

s 

H - [ h j (T) + h j + C j] Y j + I “ ' “ l»-22] 

j“l 

The equations given above are transformed to the curvilinear co- 
ordinate system depicted in Figure H-2. They are then numerically 
integrated, using Runge-Kutta integration techniques to give the 
detailed flow field solution. The details of the transformation 
and integration techniques are given in Reference 1 and will not 
be repeated here. It is worth noting, however, that convergence 
of the integration scheme (and location of the body-point along 
a ray) is determined by mass-balancing. 



2 . 


The dimensional mass flux into the shock between the body symmetry 
axis and a given X-station, M, is 


;<M) = (2 




r SH^ 


[“= 


- u 


cct 0 


SH 


dr 


[H-23] 


where u^, are components of the velocity vector, is the 

freest ream gas density, and o is the bow shock angle with respect 
to the centerline at the given r-station. The exponent, u), is 
zero for two-dimensional flow, and one for axisymmetric flow- In 
the computational procedure, the mass flux is nondimensionallzed, 
giving 

/ D (w + 1) 

*nondim (M) “ *SH (M >/ P “ SH 

where Rg^ is the bow shock radius at the centerline. Therefore, 
at a general point on a ray 


^(M) = ^ SH 



(2rrr) ^pudy 


[H-25] 


Inward integration is continued until i^(M) is less than some small 
inpi.t constant, called (M) . This implies that the flow has 

reached the body* Because all of the rays do not always reach the 
same distance from the shock wave before termination of the inte- 
gration, an irregular (e.g. wavy) body surface may result. In 
such cases, iteration of the procedure with adjustments to the 
input shock shape becomes necessary. 

The Inviscid Supersonic Calculation 

The supersonic inviscid flow field is calculated by an adaptation 
of the method of characteristics. The method used has proven to 
be convergent within reasonable computer times. 1 The derivation 
of the characteristic equations is given below. A detailed dis- 
cussion of the computer techniques applied to integrate these 
equations is given in Reference 1. 

Derivation of the Characteristic Equations: 
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For an inviscid, compressible, reacting supersonic flow the con- 
servation equations may be expressed as follows. 


Conservation of Mass 


u 3p 3u / 3<f> , U) \ 

t = - — - uH 1 + — sin 

p 3s 3s \3n r 


Conservation of Momentum 


2 ^ = _ l ap 

3s p 3n 

3u = _ 1 3P 
U 3s p 3s 

Conservation of Energy 

3h = 3u 

3s 3s 

where 

<j> = flow angle 
r = the i idial coordinate 

n = a coordinate normal to the streamline, S, of interest 
s = the streamline coordinate 

! 0, for 2-dimensional flow 
1, for axisymmetric flow 
The mixture static enthalpy is given by 
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N • K 
vib 


■ E 

j=i 
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J 3 


E 

k=l 


N + k 
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[H-27] 

[H- 28 ] 


[H-29] 


[H-30 ] 
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Here, h' denotes all of the species enthalpy except that caused 

by vibration, and the second term denotes the vibrational contri- 
butions to mixture enthalpy. The variation of concentration of 

the j tb chemical species along a streamline, s, as given by the 
two following equations : 


N t 


u K/3s) ■ rE Q u 

J-l 


j = N + 1, ..., N 
c s 


3g 1 

u ^r = w j 


[H-31] 

[H-32] 


The equation of state for the gas mixture i9 
N 


= p rt£ 


'j 


j-1 


The derivation of the characteristic equations, from this system, 
may be started by eliminating as many derivatives as possible 
from Equations [H-26] through [H-32]. Consider, first, Equations 
[H-28], [H-29], and [H-30]. Differentiation of Equation [H-30], 
with respect to s, yields: 


[H-33] 
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Noting that 


s 

E 

j=i 


y fill il . c II 

Y j 3T 3s p 3s 


gives the result that 


3h 

3s 


_ 3T , V' V.' !ll , V' / 3e i 3y N +k \ 

- C pal + E h jas t S(v?r t * k -5-) 


[H-35] 


Substitution of Equations [H— 28] and [H-29] into Equation [H-30] 
yields 


N 


1 3P = 3h = c 
p 3s 3s p 


3 Y, 


vib 


3e, 


n + y h - -i+ y Y 

3s j 3s Z-f 'N +k 3s 

j k=l c 

Now from the equation of state, e.g. Equation [H-33], one obtains 


[H-36 ] 



Now if Equation [H-26] and [H-28] are combined, the result is, 

3p _ 1 _3P _ / 3<)) , to sin j> \ 

3s u 2 " 3s p \3n r j 


[H-37] 


[H-38] 


If Equation [H-37] is substituted into Equation [H-38], the re- 
sult is 
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3P 

3s 
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( 3<t> , ui sin $ 
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[H-39 ] 


H-13 



If Equation [n**36] is solved for — , and the resulting expression 

3 3 

is substituted into Equation [H-39 ] , the result is, 




| 3P 3 j) __ u) sin 

3s 3n r 
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[H-40] 


Note that, of the original system of equations, all except Equa- 
tion [H-27] have been used in the derivation of Equation [H-40]. 
Therefore, the original set of equations governing the inviscid 
supersonic flow field has been reduced to a set of two equations, 

involving four derivatives, and 


The new subsystem Equations [H-27] and [H-40] must be supple- 
mented with two additional equations, to solve for the four un- 
knowns just given. The two supplementary equations selected by 
Curtis and Strom 1 are 

dr) + ~~ ds = d [H— 41] 

3 r t 3 s 

dn + ds = dP [H-42] 

3rj 3s 

Since the characteristic curves are defined as the loci of in- 
determinate flow derivatives, solving the set composed of Equa- 
tions [H-27], [H-40], [H— 41] , and [H-42], and locating the regions 
of indeterminacy gives the characteristic lines of the flow field. 
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This procedure can be illustrated by grouping the four equations 
together and solving for the pressure gradient. This will yield 
explicit expressions for the characteristic network in the flow 
field. Therefore, application of Cramer's Rule to this set of 
equations yields the following expression for the pressure 
gradient : 
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Where A 1 is the coefficient of — in Equation [H-40], and B 

ds 

the righthand side of Equation [H-40]. 


is 


If — is to be indeterminant , then the numerator and denominator 

3s 

must vanish separately. From the vanishing of the denominator, 
the following relation results: 


*l-± 

3s V pA 1 ^ 


tan [5 


[H-44] 


Which gives two differential equations for two sets of character- 
istic curves. Along these curves, the numerator must also vanish, 
yielding. 


pu 2y/ pA^u 2 dP ± d(j) = - 


B j ds 

-/p A*U 2 


[H-45 ] 


The compatibility equation may be obtained by using the character- 
istic angle, 3, and a length element, d£, defined as 


d ^± 


ds 

cosB 


[H-46] 


where the ± notation indicates the direction of the characteristic 
line. 
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The resultant compatibility equation is 


dP ± dq> - sin<J) + F j sing d; + [H-47 ] 

in which all nonequilibrium effects are expressed in the func- 
tional 
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Along the characteristic curves the partial differential equa- 
tions reduce to ordinary differential equations. 

The compatibility equation and the characteristics equations are 
normalized using the same reference quantities as in the subsonic 
flow field. These are R cu , p , u , p u 2 , RT , and RT /(MW) . 

bh 00 00 00 OO oo 00 00 

The details of the normalization and computational procedures are 
given in Reference 1. 


Reference 1 : 

J. T. Curtis and C. R. Strom: Computations of the Nonequilibrium 

Flow of a Viscous } Radiating Fluid About a Blunt Axisymmetric 

Body. AFFDL TR-67-40, Volume 1, June 1967. 
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APPENDIX I 


MULTILAYER INSULATION 
CONDUCTIVITY EVALUATION 


C. Webb 


February 1, 


1972 


To evaluate insulation performance accurately, two applicable data 
sources on multilayer insulation were obtained. The first source 
represented laboratory guarded hot plate test data (Ref 1), the 
second source represented thermal conductivity test data from an 
actual hardware mockup where fiberglass standoffs penetrated the 
insulation blanket and seams and joints were present (Ref 2) . 
Figure 1-1 presents a comparison between these two sets of data. 

In addition, an analytical curve fit is shown for each set where 
the expression representing the thermal conductivity considered 
both the linear conduction and the radiation associated with 
multilayer insulation. The expression, therefore, includes the 
influence of the mean insulation temperature and, in addition, 
the nonlinear influence of the actual boundary temperatures. The 
expression used for the multilayer insulation follows: 

k ■ aI M + b (V + V) ( T H + T c) [I ’ 11 

where, 


k = effective thermal conductivity 
a & b = influence coefficients 

T^ = insulation mean temperature 

T^ = hot boundary temperature 

T^ = cold boundary temperature 

and 

( t h - T c) 

where 


Q 


leak 


A 


t 


Blanket heat leak 
= insulation surface area 
= insulation thickness. 


U-2] 


The importance of this expression is that it more precisely de- 
termines the thermal conductivity as a function of boundary tem- 
perature and thus insulation thickness variations. 

For the turbopause mission thermal analyses, the Skylab conductivity 
data was selected to determine the baseline influence coefficients 
of Equation [1-1] and thus the multilayer insulation performance. 
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Figure 1-1 Multilayer Insulation Thermal ConduatiiKty 
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APPENDIX J 


THERMAL SENSITIVITY ANALYSIS 


C. Webb 


February 1, ]972 


Perhaps the biggest uncertainty in determining the amount of 
internal heat required for thermal control was the thermal con- 
ductivity of the multilayer insulation blanket. Since the methods 
of installation and fabrication can have considerable effect on 
the efficiency of the insulation, a sensitivity analysis was con- 
ducted to determine to what extent lack of knowledge of the actual 
blanket conductivity could affect the heater power requirement. 
These results are presented in Tigure J-l. It is expected that 
effective insulation conductivity will fall somewhere between 
0.0002 and 0.002 W/m - °K. The response of heater power to this 
order of magnitude variation is approximately a factor of 2 (from 
7 to 16 watts). This means that the final decision on the number 
of radioisotope heaters required for the probe should be made on 
the basis of a thermal vacuum test. This test would be relatively 
simple and could be performed on the actual flight article, thus 
precluding the expense of fabricating a thermal test prototype 
and establishing maximum confidence in thermal performance. Add- 
itional uncertainties in the thermal analysis were the. optical 
< properties of the beryllium heat shield, i.e., solar absorptivity 
and emittance. Since parts of the shield will be plated with 
platinum, a sample platinum-plated beryllium chip was prepared 
and laboratory optical property measurements performed. The 
platinum plating was accomplished by a thermally decomposable 
paint which, after heat treatment, leaves a platinum coating. A 
laboratory sample of cross rolled beryllium sheet was supplied 
and prepared by Englehard Minerals and Chemicals, Newark, New 
Jersey. Emissivity measurements performed by Martin Marietta 
indicated a value of e s 0.25 for platinum-coated beryllium. 

These measured values are not expected to be the definitive 
answer for probe exterior emissivity. Considerable differences 
can occur, depending on the degree of polish and plating char- 
acteristics . 

Analytical results showing the influence of the probe shell optical 
properties on the internal heater power requirement are presented 
in Figure J-2. These results show that the heater power is not 
particularly sensitive to variations in shell emissivity for 
values 0.20 or greater. Lower emissivity values, however, could 
save significant heater power but coincidently enlarge the uncer- 
tainty associated with the design due to the extreme sensitivity 
of heater power to variations in the probe shell emissivity in 
this range. For this reason, it is recommended that some effort 
be made, primarily through control of surface finish (polish) to 
ensure that the probe shell actually has an emissivity of 0.2 or 
greater . 
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APPENDIX K 


GROSS MISSION ANALYSIS FOR EARTH - JUPITER 
MISSIONS AND A JUPITER PROBE FOR 1976-1983 


L. Anderson 
D. Vogt 


September 28, 1971 


A parametric analysis has been completed to aid in the selection 
of Earth Launch Date (LD) , Jupiter Arrival Date (AD) and inertial 
flight path angle of probe at entry (y ) for the 1976, 1977, 1978, 

1979, 1980, 1981-1982, and 1983 launch windows. This analysis is an 
addendum to the data presented in Ref 1. 

The following flight parameters were analyzed over the complete 

range of LD (1976 through 1983), AD and y (7.5 to 67.5 deg): 

E 

Launch energy, C3 (km/sec) 2 

Equatorial declination of departure v , DLA (deg) 

HP 


Sun-Ear th-vehicle angle at Jupiter arrival, SEV (deg) 

Navigation angle at Jupiter arrival, ZAE (deg) 

Sun-Jupiter - v angle at Jupiter arrival, ZAP (deg) 
HP 


Probe illumination angle at entry, X (deg) 

Probe-Earth communication angle at entry, (deg) 

Minimum entry angle consistent with 20° mask angle, y . (deg) 

min 

These parameters are discussed in detail in Ref 1 and 2. 

The constraints on these parameters follow. 

Soft Constraint s Caused by 

C3 < 130 Launch vehicle capability (T3D-Centaur ) 

( DLA j > 2 Navigation for first midcourse correction 

| DLA | < 36 Range safety 

Doppler navigation approaching Jupiter 


ZAE + 90 
ZAP <150 
X < 70 

hi < 20 


Minimum p^obe y requirements 
E 

Minimum probe lighting condition 

Favorable probe - Spacecraft-Earth 
communication geometry 
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Hard Constraints Caused by 


C3 < 170 


Launch vehicle capability (T3D-7 Centaur) 


DLAj < 40 Range safety 

SEV| > 15 Sun obstruction at Jupiter arrival 


A. MISSION ANALYSIS 


For all launch windows, there is a tradeoff that must be made 

between probe y and launch window duration. If science, aero- 
L 

heating, or communications require that the probe enter with 
lower y , such as less than 30 degrees, and still maintain light- 
ing conditions, then the window of allowable ADs is restricted to 
later dates, reducing the available launch window as illustrated 

in Figure K-l. This relationship between minimum possible y , AD, 

L 

and launch window is illustrated in the C3 plots (Figure K-2) . 
Reference 3 states that the minimum launch window requirement for 
the Titan IIID/5 segment vehicle is 20 days, and the minimum for 
the Titan IIID/7 segment vehicle is 30 days. These requirements 
are indicated on Figure K-l. 

There is also a tradeoff to be made between y r , and time of flight, 

since smaller probe y r angles require later Jupiter ADs. This 

£> 

tradeoff is also approximately illustrated on Figure K-l. For all 
launch windows, the time of flight can vary from about 450 to 1250 
days . 

The only other AD limitation common to all AD windows is the sun 
obstruction constraint. The sun obstruction occurs for two 40-day 
intervals during the 750-day AD window, but is easily avoided by 
small changes in Earth-Jupiter targeting. 

The only significant differences in the seven LD windows considered 
are the presence of range safety constraints in 1977, 1978, and 
1979 and a midcourse navigation constraint in 1981-1982. Assuming 
that all of these constraints are hard constraints, the launch 
windows fall into the order of most favorable to least favorable 
as follows: 
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LD Window 

Minimum y ^ Possible* 

Comments 

1976 

0 

No launch constraints 

1983 

0 

No launch constraints 

1980 

0 

No launch constraints, narrower 
LD window 

1977 

16.5 

Partly constrained by range 
safety 

1979 

20 

Partly constrained by range 
safety 

1978 

26.5 

Partly constrained by range 
safety 

1981-82 

0 

Completely constrained by 


midcourse navigation 

^Approximate , assuming a 20-day launch window requirement and a 
launch vehicle capability of C3 not greater than 130. 


PRESENTATION OF DATA 


Figures K-2 through K-15 present the new data generated for ^his 
analysis. This form of data presentation is made possible by 
recent modifications tu the computer programs GE0M1 (written by 
A. Barger) and C0NT0UR (written by F. Peterson). Additional 
parametric pres ntations of this form can be generated as needed. 

The C3 contour graphs are similar to these prese ited in Ref 1, 
with the following exceptions. The basic data and hardcopy plots 
are generated by single computer run instead of generating 
rough data first and then having graphs hand plotted. The graphs 
include the parameters ZAE and ZAP which have not been previously 
plotted. 

The X, \p plots are similar to data presented in Ref 1, but cover 
ot^er LD windows in addition to 1978. Again the plots are machine- 
generated instead of hand-drawn. The data on each X, 4 , plot cor- 
responds to a single launch date at the center of the LD window 
as noted in the X-axis label. TV 0 X, \p plots, Figures *.<-5 and 
K-6, were generated for whe ’977 window, correspon 'ir z o two 
different launch dates. A U-day change in LD produced only small 
changes (less than one degre^) in the X and ip contours. 
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CONCLUSIONS 


Considering the 1978, 1979, and 1930 launch windows, the 1978 is 
least attractive for flight path angles less than 30 degrees. 

This is due to the range safety constraint that eliminates a 
significant portion of the launch window. For example, if the 
probe is targeted to a y of 20 degrees, then the C3 requirements 

are 150 (km/sec) 2 for a 20-day launch window and 200 (km/sec) 2 
for a 30-day launch window. Either 1979 or 1980 presents a 
better alternative, if the probe must enter at low y £ . 
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CONCLUSIONS 


Considering the 1978, 1979, and 1980 launch windows, the 1978 is 

least attractive for flight path angles less than 30 degrees. 

This is due to the range safety constraint that eliminates a 

significant portion of the launch window. For example, if the 

probe is targeted to a y of 20 degrees, then the C3 requirements 

E 

are 150 (km/sec) 2 for a 20-day launch window and 200 (km/sec) 2 

for a 30-day launch window. Either 1979 or 1980 presents a 

better alternative, if the probe must enter at low y . 
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Fi<?. F-J Entry Angle, Launch Window and Flight Time Tradeoff, ±976 to 1983 
















Arrival Date 


8/ 17/Ml 
7/18/81 
8/18/81 
5/19/81 
4/19/81 
3/20/81 
2/18/81 
1/19/81 
12/20/80 
11/20/80 
10/21/80 
9/21/80 
8/22/80 
7/23/80 
6/23/80 
5/24/80 
4/2<*/8Q 
3/25/80 
2/24/80 
1/25/80 
12/26/79 
11/26/79 
10/27/79 
9/27/79 
8/28/79 
7/29/79 
6/29/79 
5/30/79 
4/30/79 
3/31/79 
3/1/79 
1/30/79 
12/31/78 
12/1/78 
11/1/78 


• 90 


ZAP ■ 90, r m - 0 



) I SEV | < 15° 
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Fig. K-9 Jupiter Turbopause Probe Study 
1979 Launch, Type I and II 
(TF at X, Y * (0,0) is 454 Days) 
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Fig. K-10 Jupiter Turbopause Probe Study 
1979 Launch, Type I and II 
(TF at X, X = (0,0) is 454 Days ) 
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Fig. X-J2 Jupiter Turbopause Probe Study 
1980 Launch, Type I and II 
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Fig. K-13 Jupiter Turbopause Probe Study 
1981-82 Launch , Type I and II 
(TF at X, Y = (0,0) Is 451 Days ) 


K-17 



Arrival Date 


12/15/85 
11/1* '85 
10/1 
9/16/85 
8/17/85 
7/18/85 
6/18/85 
5/19/85 
4/19/85 
3/20/85 
2/18/85 
1/19/85 
12/20/84 
11/20/84 
10/21/84 
9/21/84 
8/22/84 
7/23/84 
6/23/84 
5/24/84 
4/24/84 
3/25/84 
2/24/84 
1/25/84 
12/26/83 
11/26/83 
10/27/83 
9/27/83 
8/28/83 
7/29/83 
6/79/83 
5/30/83 
4/30/83 
3/31/83 
3/1/83 



Gama (LD ■ 1/14/82) 


Fig. F-74 Jupiter Turbopauae Probe Study 
1981-82 Launch, Type I and II 
(TF at X, I - (0,0) la 461 Days) 
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